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Abstract 

This talk at Fradkin conference is devoted to brief description of 
latest results in two topics I worked on during last years. 

The multiresolution analysis and fast wavelet transform became a 
standard procedure for pattern recognition and I describe how it has 
been used for analysis of high energy nucleus-nucleus interactions^. 

Latest QCD results on multiplicities in quark and gluon jets are 
discussed and confronted to experimental dataP]. 

1 Introduction 

This talk is given at Fradkin memorial conference on Friday, 9 June 2000. 
Friday was always a traditional day of seminars on quantum field theory 
in our department, so-called Fradkin seminars. Two talks on wavelets were 
given at this seminar last year just before Fradkin was put in a hospital. 
Therefore I decided to talk here about wavelet application in particle physics 
and gave to organizers of the conference the title containing the first part of 
the above one. However soon I learned that the topic on QCD (the second 
part of the present title) would fit the conference schedule better and asked 
organizers for replacement. It was shifted to QCD session but with an old 
title on wavelets. Thus I decided to give two talks at one, and it explains 
how these two topics appear here together. 



2 Wavelets 

First I learned about wavelets from Pete Carruthers in 1993. He applied 
them m, H, |H for analysis of some scaling cascade models used, in particular, 
in multiparticle production modelling. It was briefly described in our review 
paper Then I proposed to use wavelets for pattern recognition in high 

1 Other applications can be found, e.g., in Web site www.awavelet.com. 
2 Detailed review of this subject is given in |Q 
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energy nucleus-nucleus collisions, and it was applied to experimental data 

Wavelets became a powerful mathematical tool in many investigations. 
They are used in those cases when the result of the analysis of a particular 
signal^ should contain not only the list of its typical frequencies (scales) but 
also the knowledge of the definite local coordinates where these properties 
are important. Wavelets form a complete orthonormalized system of func- 
tions with a finite support by using dilations and translations. That is why 
by changing a scale (dilations) they can describe the local characteristics of 
a signal, and by translations they cover the whole region in which it is stud- 
ied. Due to the completeness of the system, they also allow for the inverse 
transformation to be done. The locality property of wavelets leads to their 
substantial superiority over Fourier transform which provides us only with 
the knowledge of global frequencies (scales) of the object under investigation 
because the system of functions used (sine, cosine) is defined on the infinite 
interval. 

High energy collisions of elementary particles result in production of many 
new particles in a single event. Each newly created particle is depicted kine- 
matically by its momentum vector i.e. by a dot in the three-dimensional 
phase space. Different patterns formed by these dots in the phase space 
would correspond to different dynamics. To understand this dynamics is a 
main goal of all studies done at accelerators and in cosmic rays. Especially 
intriguing is a problem of the quark-gluon plasma, the state of matter with 
deconfined quarks and gluons which could exist during an extremely short 
intervals of time. One hopes to create it in collisions of high energy nuclei. 
Nowadays, the data about Pb-Pb collisions are available where in a single 
event more than 1000 particles are produced. We are waiting for RHIC ac- 
celerator in Brookhaven and LHC in CERN to provide events with up to 
20000 new particles created. Therefore the problem of phase space pattern 
recognition in an event-by-event analysis becomes meaningful. It is believed 
that the detailed characterization of each collision event could reveal the rare 
new phenomena, and it will be statistically reliable due to a large number of 
particles produced in a single event. 

When individual events are imaged visually, the human eye has a ten- 

3 The notion of a signal is used here for any recorded information about some processes, 
objects, functions etc. 
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dency to observe different kinds of intricate patterns with dense clusters 
(spikes) and rarefied voids. However, the observed effects are often dominated 
by statistical fluctuations. The method of factorial moments was proposed 
[||] to remove them but it is hard to use in event-by-event approach. The 
wavelet analysis avoids smooth polynomial trends typical for the statistical 
component. It was first applied || to analyze the individual high multiplic- 
ity event of Pb-Pb interaction at energy 158 GeV per nucleon. With emul- 
sion technique used in experiment the angles of particle emission are often 
measured only, and the two-dimensional phase space is considered therefore. 
The experimental statistics is rather low but acceptance is high and homo- 
geneous that is important for proper pattern recognition. To simplify the 
analysis , the two-dimensional target diagram representing the polar and az- 
imuthal angles of created charged particles was split into 24 one-dimensional 
functions representing the polar angle distribution of these particles in 24 
azimuthal angle sectors of tt/12 and in each of them particles were projected 
onto the polar angle 9 axis. Thus one-dimensional functions of the rapidity 
distribution of these particles in 24 sectors were obtained. Then the wavelet 
coefficients were calculated in all of them and tied up together (continuous 
MHAT wavelet was used). The resulting pattern showed that many parti- 
cles are concentrated close to some value of the polar angle i.e. reveal the 
ring-like structure in the target diagram. The interest to such patterns is 
related to the fact that they can result from the so-called gluon Cherenkov 
radiation [[?], || or, more generally, from the gluon bremsstrahlung at a finite 
length within a quark-gluon medium (plasma, in particular). More elabo- 
rated two-dimensional analysis was done recently |J and confirmed these 
conclusions with jet regions tending to lie on some ring-like formations. The 
jet-like substructure of the event becomes more pronounced, and ring-like 
correlations of jettty regions are noticeable. With higher statistics, one can 
learn if the angular distribution of these rings corresponds to theoretical ex- 
pectations. It is due to wavelet analysis that for the first time the fluctuation 
structure of an event is shown in a way similar to the target diagram on the 
two-dimensional plot. 



Previously, some attempts [10, 11, 12 to consider such events with differ- 



ent methods of treating the traditional projection and correlation measures 
revealed just that such substructures lead to spikes in the angular (pseudora- 
pidity) distribution and are somewhat jetty. Various Monte Carlo simulations 
of the process were compared to the data and failed to describe this jettiness 
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in its full strength. More careful analysis |13], [13] of large statistics data 
on hadron-hadron interactions (unfortunately, however, for rather low mul- 
tiplicity) with dense groups of particles separated showed some "anomaly" 
in the angular distribution of these groups awaited from the theoretical side. 
Further analysis using the results of wavelet transform are needed to check 
this conclusion in high multiplicity nucleus-nucleus interactions when many 
events of this kind become available. 



3 QCD 

Here I briefly describe recent advances in theoretical understanding of multi- 
plicity distributions of quark and gluon jets. The extended survey with more 
detailed comparison to experimental data can be found in the recent review 
paper JT5|| . 

The progress in experimental studies of properties of quark and gluon 
jets is very impressive. Therefore the study of the energy evolution of such 
parameters of multiplicity distributions of jets as their average multiplicities 
and widths becomes possible. It is well known that the average multiplicities 
of quark and gluon jets increase quite fast with energy but their ratio has a 
much slower dependence. 

The perturbative QCD provides quite definite predictions which can be 
confronted to experiment. In brief, the results can be summarized by saying 
that the energy dependence of the mean jet multiplicity can be perfectly fitted 
but the ratio of gluon to quark jet multiplicities can be described within the 
precision of 15-20% only. Moreover, one can understand why next-to-leading 
approximation is good enough for describing the energy dependence, but it 
is not quite satisfactory yet for the ratio value. I show this by presenting the 
analytical expressions. For the corresponding Figures, I refer the reader to 
the review paper [T5| . 



The theoretical asymptotical value of the ratio of average multiplicities 
equal 2.25 is much higher than its experimental values, which are in the 
range from 1.05 at comparatively low energies of T resonance to 1.5 at Z° 
resonance. The next-to leading order (NLO) corrections reduce this ratio 
from its asymptotical value by about 10% at Z° energy. The NNLO and 
3NLO terms diminish it further and show the tendency to approximate the 
data with better accuracy. The computer solution of QCD equations for 
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the generating functions has shown even better agreement with experiment 
not only on this ratio but on higher moments of multiplicity distributions as 
well. Being perfect at Z° energy, the agreement in the ratio is not as good 
at lower energies where the theoretical curve is still about 15-20% above the 
experimental one. In other words, the theoretically predicted slope of the 
ratio of multiplicities in gluon and quark jets is noticeably smaller than its 
experimental value. Nevertheless, one can speak about the steady conver- 
gence of theory and experiment with subsequent improvements being done. 
Moreover, it is even surprising that any agreement is achieved in view of the 
expansion parameter being extremely large (about 0.5) at present energies. 

The importance of studying the slopes stems from the fact that some of 
them are extremely sensitive (while others are not) to higher order perturba- 
tive corrections and to non-perturbative terms in the available energy region. 
Thus they provide us with a good chance to learn more about the structure 
of the perturbation series from experiment. 

In the perturbative QCD, the general approach to studying the multiplic- 
ity distributions is formulated in the framework of equations for generating 
functions. Therefrom, one can get equations for average multiplicities and, in 
general, for any moment of the multiplicity distributions In particular, 



two equations for average multiplicities of gluon and quark jets are written 



as 



(n G (y)) = J dxil[K%{x){(n G {y + lux)) + (n G (y + ln(l - x)) - (n G (y))) 

+n f K^(x)((n F (y + lux)) + (n F (y + ln(l - x)) - (n G (y)))}, (1) 

(n F (y))' = [ dx^K^(x)((n G (y + \nx)) + (n F (y + \n(l-x))-(n F (y))). (2) 



Herefrom one can learn about the energy evolution of the ratio of multiplic- 
ities r and of the QCD anomalous dimension 7 (the slope of the logarithm 
of average multiplicity in a gluon jet) defined as 

fog) K?)' n / v v 

{op) ' 7 = Wi = iHnG)) ■ (3) 

Here, prime denotes the derivative over the evolution parameter y = hx(pQ/Qo), 
p, are the momentum and the initial angular spread of the jet, related to 
the parton virtuality Q = pQ/2, Qo=const, K's are the well known split- 
ting functions, (n G ) and (n F ) are the average multiplicities in gluon and 
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quark jets, (no)' is the slope of (%), rif is the number of active flavours. 
The perturbative expansion of 7 and r is written as 

7 = 70(1 - 0170 - a 2 7o - as7o) + °(To)i ( 4 ) 
r = r (l - ri7o - r 2 7o - r 3 ^) + 0(7q), (5) 



where 70 = \j2N c as/iT, as is the strong coupling constant 

2tt 



as = 75— 



1 Aln(2y) 



0(|/- 3 ), (6) 



A) = (HiVc - 2n/)/3, A = (51iV c - 19n/)/3,r = iV c /C F , and in QCD 
A^ c = 3 is the number of colours, Cp = 4/3. 

The limits of integration in eqs. ([[]), (@) used to be chosen equal either 
to and 1 or to e~ y and 1 — e~ y . This difference, being negligibly small at 
high energies y, is quite important at low energies. Moreover, it is of physics 
significance. With limits equal to e~ y and 1 — e~ y , the partonic cascade 
terminates at the perturbative level Qo as is seen from the arguments of 
multiplicities in the integrals. With limits equal to and 1, one extends the 
cascade into the non-perturbative region with low virtualities Qi ~ xpQ/2 
and Q 2 ~ (1 — x)pQ/2 less than Qo/2. This region contributes terms of 
the order of e~ y , power-suppressed in energy. It is not clear whether the 
equations and LPHD hypothesis are valid down to some Qo only or the non- 
perturbative region can be included as well. 

Nevertheless, the purely perturbative expansion (^), (|5p with constant 
coefficients a^, and energy-dependent 70 is at work just in the case of limits 
and 1. The values of a^, Tj for different number of active flavors rif are tabu- 
lated in [0 . At Z°-energy the subsequent terms in (||]) diminish the value of 
r compared with its asymptotics r = 2.25 approximately by 10%, 13%, 1% 
for nf = 4 getting closer to experiment. However the theoretical value of r 
still exceeds its experimental values by 15-20%. 



The energy dependence of mean multiplicities can be obtained |T8| [17 
from the definition (^) by inserting there the value of 7 (|]) and integrating 
over y. Keeping the terms as small as y^ 1 at large y in the exponent, one 
gets [{r7| the following expressions for energy dependence of multiplicities of 
gluon (G) and quark (F) jets 



(n G > = Ky- a ^ exp [2C^ + 5 G {y)\ , (7) 
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with K an overall normalization constant, C = J4N c /(3q } and 



C 



2a 2 C 2 + §[ln(2 2/ ) + 2] 



+ - 



C 2 



«3C 2 -^#N2 2/ ) + l] 



(n F ) = —y 
r 



ai c 2 



exp \2C^/y + #f 



with 



C C 2 r 2 
S F {y) = S G {y) + — n + — (r 2 + 



(8) 



(9) 



(10) 



It happens that 2NLO and 3NLO terms (contributing to y~ 1 ^ 2 and y~ l terms 
in the exponent) are almost constant at present energies and do not change 
the energy dependence prescribed in NLO approximation. It explains why 
the energy dependence is well fitted by both NLO and 3NLO formulas while 
2NLO correction to the ratio r is large and important. 

The rather small difference in r values results in quite noticeable dis- 
agreement of the slopes r'. Theoretical estimates can be shown ]17| to be 
quite predictive for the ratio of the slopes of multiplicities but it is much less 
reliable to use the perturbative estimates even at Z°-energy for such quan- 
tities as the slope of r or the ratio of slopes of logarithms of multiplicities 
(the logarithmic slopes). Much higher energies are needed to do that. Thus 
the values of r and/or of the logarithmic slopes can be used to verify the 
structure of the perturbative expansion. 

I demonstrate it here on the example of the slope value. The slope r' 
is extremely sensitive to higher order perturbative corrections. The role of 
higher order corrections is increased here compared with r because each nth 
order term proportional to 7q gets an additional factor n in front of it when 
differentiated, the main constant term disappears and the large ratio r 2 /ri 
becomes crucial: 



2ro7n 

n 



3r 3 

n 



+ fi 1 ) 7 o + 0(7o 3 ) 



(11) 



where the relation 7g w — £>7q(1+.Bi7q) has been used with B = (3 /8N c ; B\ = 
/3i/4:N c /3o. The factor in front of the bracket is very small already at present 
energies: Br^ri pd 0.156 and 70 ~ 0.5. However, the numerical estimate of r 
is still indefinite due to the expression inside the brackets. Let us note that 
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each differentiation leads to a factor as or 7g, i.e., to terms of higher order. 
For values of 7*1, r2, r% tabulated above (n/ = 4) one estimates 2r2/r\ m 4.9, 
(Srs/rij+Bi ~ 1.5). The simplest correction proportional to 70 is more than 
twice larger 1 at energies studied and the next one is about 0.4. Therefore 
the ever higher order terms should be calculated for the perturbative values 
of r' to be trusted. The slope r' is equal to for a fixed coupling constant. 

The higher order terms are also important for the moments of the multi- 
plicity distributions [19| . The normalized second factorial moment F2 defines 



F^as — q! ^2,as ~ A - (12) 



the width of the multiplicity distribution. 

The asymptotical (70 — > 0) values of F 2 G and F 2 F are different: 

4 ci 7 
- F = - 

g ) 2, as j ' 

At Z° energy the widths of the distributions are smaller 

F 2 G « 1.12, F 2 F « 1.34. (13) 

but still much larger than their experimental values 1.02 and 1.08, corre- 
spondingly. The rather large difference of the perturbative ( |13"D and exper- 
imental values at Z° indicates that moments of the distributions should be 
sensitive to corrections. The conclusions about the third moments are sim- 
ilar. Nonetheless, the computer solution of QCD equation happened to be 
quite successful in fitting experimental data even for higher moments and 
their ratios H q introduced in p0| . It shows that the role of conservation laws 
treated approximately in the analytical approach and accurately accounted 
in computer calculations becomes more important for higher moments. 

Thus it is shown that the analytical approach is quite successful in demon- 
strating that all features of QCD predictions about multiplicities of quark 
and gluon jets correspond to the general trends of experimental data. Some 
disagreement at the level of 15-20% is understandable due to incomplete ac- 
count for the energy-momentum conservation in such an approach. Further 
accurate computer solutions are needed to check if these trends persist at the 
higher precision level. 
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